复合推进剂侵蚀燃烧特性的κ-ε参数表征

本文从推进剂侵蚀燃烧机理分析出发,综合片型装药发动机中止燃烧试验与相应的燃烧室内流场N-S方程组STMPLE方法数值分析的特点,导出一种由燃面上k~2/ε和p参数表述推行剂侵蚀特性的新方法,并成功地应用于4B丁羟复合药侵蚀特性的研究中.     

嵌金属丝推进剂燃烧起始阶段增速特性研究

考虑了使用嵌金属丝装药的发动机工作过程中的多种因素,建立了金属丝与药柱、燃气间相互作用的控制方程组,并对方程组进行了数值求解,得到了嵌金属丝推进剂在燃烧起始阶段,增速比随时间的详细变化过程。结果表明,金属丝的热扩散率、熔点、直径和黑度以各自不同的方式影响金属丝的增速效果,熔点相对较低的金属丝可能会在内弹道曲线的初始阶段形成一个轻微的压强峰,固有燃速较低、着火温度较低的推进剂对嵌入金属丝较为敏感。文中数值方法计算量小,结果可靠,适于工程应用。     

高旋转数内冷通道换热实验技术及验证

通过内冷通道换热理论分析及实验系统优化,获得了使内冷通道内气体压力保持在500kPa以上的高旋转数实验技术.旋转数和雷诺数的范围分别为0~2.08,104~7×104,极大满足了旋转数与雷诺数同时覆盖真实发动机涡轮叶片工作参数的实验需求.证实了雷诺数和旋转数对内冷通道换热的影响可以解耦.而应用该实验技术进行同等旋转数实验时,内冷通道热损失占总加热量的比例和实验误差分别降低到18%和±10%.     

THE PREDICTIONS OF CONVECTIVE HEAT TRANSFER ON TURBINE BLADE AIRFOIL BY USING LOW－REYNOLDS NUMBER TURBULENCE MODEL

ＴＨＥＰＲＥＤＩＣＴＩＯＮＳＯＦＣＯＮＶＥＣＴＩＶＥＨＥＡＴＴＲＡＮＳＦＥＲＯＮＴＵＲＢＩＮＥＢＬＡＤＥＡＩＲＦＯＩＬＢＹＵＳＩＮＧＬＯＷ－ＲＥＹＮＯＬＤＳＮＵＭＢＥＲＴＵＲＢＵＬＥＮＣＥＭＯＤＥＬＴＨＥＰＲＥＤＩＣＴＩＯＮＳＯＦＣＯＮＶＥＣＴＩＶＥ...     

点阵单元类型和排布角度对传热特性的影响

基于瞬态液晶测试技术和数值仿真,在相同孔隙率条件下,研究了桁架型点阵结构的单元类型和排布角度对矩形通道壁面对流传热特性的影响.结果表明:点阵单元类型对壁面传热性能有决定性的影响,其中交叉排布的体心立方(body centered cubic-0,BCC-0)和Kagome-0点阵归因于出色的传热强化效果,雷诺数在5 7...     

光敏推进剂激光烧蚀过程建模与仿真

相较于传统大卫星,微小卫星具有结构紧凑、质量轻便和成本低廉的特点。然而,受功率和质量负载的限制,微小卫星一般不装备推进系统,其航线也局限于近地轨道。为扩展微小卫星的功能,满足日益复杂的任务需求,需给其配备合适的微推进系统。固体推进系统具有结构简单、寿命长、可靠性高的优点,但无法重复启动。为得到可重复启动的固体微推进系统,设计了一种非自持燃烧的光敏推进剂,采用激光控制其燃烧。在背压为大气压的环境下,利用高速摄像机拍摄燃烧过程并记录燃速。之后,对光敏推进剂的激光烧蚀过程进行建模。分析结果表明:激光可控制光敏推进剂的燃烧,燃速与激光强度成线性关系;该光敏推进剂的最小激光点火强度为0.28 W/mm~2;燃速计算值与实测值的误差在10%以内,证明该数学模型具备工程应用价值。     

篦齿封严风阻温升特性研究

篦齿封严风阻温升效应引起的热负荷对航空发动机涡轮叶片冷气系统有着重要的影响。采用理论分析、数值计算与实验相结合的方法系统地研究了篦齿封严的风阻温升特性。首先,对篦齿封严风阻温升特性进行了理论分析,设计搭建了篦齿封严风阻温升特性实验台,建立了基于RNG（Re-Normalization Group） k-ε湍流方程的篦齿封严风阻温升数值求解模型。然后,研究了篦齿封严流场特性、泄漏特性和风阻温升特性,并将理论计算、数值仿真与实验测试结果相互对比分析,研究了压比、转速等因素对篦齿封严风阻温升特性的影响规律,揭示了篦齿封严的风阻温升效应产生的机理。结果表明：高低齿篦齿封严结构减弱了篦齿封严的透气效应,增强了篦齿封严的动能耗散,有利于降低篦齿封严的泄漏量;在所研究的工况下,转速低于2 000 r/min时,风阻温升效应较小,转速在2 000~6 000 r/min时,风阻温升随转速的升高而增大,温升值最高可达12.87 K;压比的增大会加强气流的对流换热,转速为6 000 r/min时,压比从1.1增加到1.3,温升值下降了7 K左右;风阻温升产生的主要原因是流经封严间隙的黏性气流与高速旋转的转子相互摩擦产生热量,气流吸收这部分摩擦热导致温度升高,转子转速越高,风阻温升效应越强。所研究的篦齿封严风阻温升特性为航空发动机内通道气流热负荷分析提供了理论依据。     

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.      

低压下镁铝富燃料固体推进剂燃烧性能研究

通过调节氧化剂含量、粒度级配,或加入KP、用硼粉替换部分铗铝及改变催化剂含量,研究了镁铝富燃料固体推进剂燃速和压强指数的变化规律。研究结果表明,燃速随着AP、KP粒度的减小而增加;随AP含量的增加而增加;随催化剂含量的增加而增加;随KP含量的降低而增加。压强指数随AP粒度减小呈先升高后降低再升高的趋势;粒度不同的配方随AP含量增加,压强指数变化趋势不同;催化剂含量对压强指数的影响规律也同AP粒度有关;KP的粒度变化对压强指数几乎无影响。硼粉替换部分镁铝对燃速和压强指数的影响规律与氧化剂的粒度有关。     

基于非结构网格的燃面推进算法

根据惠更斯波传播原理,构建了一种在通用CFD软件的非结构网格系统上直接计算燃面推进的数值方法。该方法将燃面对应的计算网格面作为一系列平面波源,利用网格几何关系对边界节点奇异性进行判定,通过所有平面波源影响区域的叠加求解新的燃面位置。数值实现过程中,根据节点奇异性和波源影响区域以三维矢量运算得到节点的推进矢量。结合网格的重分,在通用CFD软件环境中实现了燃面的推进,所用数值算法具有良好的精度。另外,对三维动态内弹道模拟中燃面计算的稳定性和精度等问题进行了研究。